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Abstract 



This article deals with non-adiabatic processes (i.e. processes excluded by the 
adiabatic theorem) from the geometrical (group-theoretical) point of view. An ap- 
\ proximated formula for the probabilities of the non-adiabatic transitions is derived 

in the adiabatic regime for the case when the parameter-dependent Hamiltonian rep- 
resents a smooth curve in the Lie algebra and the quantal dynamics is determined 
by the corresponding Lie group evolution operator. We treat the spin precession in 
P3 ■ a time-dependent magnetic field and the over-barrier reflection problem in a uni- 

^ . form way using the first-order dynamical equations on SU{2) and SU{1.1) group 

qh| manifolds correspondingly. A comparison with analytic solutions for simple solvable 

models is provided. 

^- 

1 Introduction. 

It is well known that probabilities of transitions induced by a time dependence of the 
Hamiltonian are suppressed if the dependence is slow. This statement formulated and 
proved by Born and Fock0 is presented in standard textbooks (e.g. by Messiah0) and 
known as the adiabatic theorem. The adiabatic theorem does not prove that the transitions 
are forbidden, it means just that the probabilities are suppressed exponentially and vanish 
to any finite order of the standard perturbation theory. Transitions of this type should 
not be discarded, however, if they result in special phenomena, even though relatively 
rare ones. The over-barrier reflection and the spin-flip in a time-dependent magnetic 
field slowly deviating from the original direction may serve as simple examples of such 
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phenomena. Usually, the adiabatic character of the process suggests a way to evaluate its 
probability, like in the quasi-classical approximation. 

Corrections to the adiabatic theorem were considered in a number of works, especially 
with application to physical systems with two non-degenerate energy levels. Dykhne]^ 
considered a Hamiltonian H{t) with two eigenstates \I'2(t) which is analytic in 

time. He found the transition probability pi2 for a system, prepared at time —oo in 
the eigenstate oo) to pass to the eigenstate \E'2 as t runs from — oo to +00. This 
probability is expressed by the formula: 

P12 ~ exp(-2r|Im j^*'(E2(t) - Eiit))dt\). (1) 

Here tc is a point in the complex time plane in which E2{t) and Ei{t) cross, and T is a 
time-scale parameter (large in the adiabatic limit) over which H{t) changes essentially. 
A rigorous derivation of Dykhne's result was given by Davis and Pechukas[^]. Suominen, 



Garraway and Stenholm|12],|13| have applied the Dykhne, Davis and Pechukas approach 



to two-level solvable models. In particular, an adiabatic behavior of the Landau and Zener 



Tj], [|l^] model was considered. It appears that for this model the Dykhne, Davis and 
Pechukas method gives the exact answer. 

From the works of Berry 0, Joye, Kunz and Pfister]^], Jaksic and Segert0,0 it be- 
comes clear that similar to Berry adiabatic phase [|l|, the transition probabilities induced 
by a time dependence of the Hamiltonian are connected with the geometry of the param- 
eter space. In particular, when the evolution of a system is described by the Hamiltonian 
of the form 

His)=n{s)-a, |n(s)| = 1, (2) 

(where n(s) is a parameter dependent unit vector and cr = (cti, (T2, cts) are Pauli matrices) 
in the first-order adiabatic perturbation theory the transition (spin-flip) W± is determined 
by the Fourier transform 0, ||^, 0, p: 

r+00 

W±r^\ exp(-2iTs)x(s)rfs|2 (3) 



X{s)=t/2\n'{s)\exp{-i^{s)), (4) 

where |n'(s)| is related with the Riemannian length element dl{s) = \n'[s)\ds of the unit 
sphere. The function q{s) is given by the integral 



?(s) = / Kg{s)ds, (5) 
Jo 

and Hg{s) is the geodesic curvature of a path n(s): 

^ n"(^)- (n;(^) xnM) . 

We note that the Hamiltonian (|^) defines a curve in the Lie algebra su(2). The cor- 
responding evolution operator belongs to the fundamental (2 x 2) representation of the 
group SU (2) and the transition probability appears to be completelly determined by the 
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geometric properties of the SU{2) group homogeneous space S'^ = SU{2)/U{1) (which 
serves as the parameter space for this particular situation). 

It is the main purpose of this work to estabhsh a relation between the probabilities of 
the non-adiabatic transitions and the geometry of the parameter space in a more general 
than the mentioned above case. Namely, we consider a situation when the commutator 
algebra of the Hamiltonian operators H{s) at different values of parameter s is closed 
for all s, constituting an arbitrary Lie algebra Q. This condition makes it possible to 
reformulate the original problem about the evolution of a quantum state in terms of a 
first-order dynamic equation on the group manifold (section 2). The adiabatic solution of 
this equation is constructed in section 3, and the integral expression for the Lie algebra 
element, determining transition probabilities, is derived. In section 4 our approach is 
applied to the spin precession in a variable external field. The adiabatic approximation 
for the spin precession 0, [Q, [|], is reconstructed. Applied to the ordinary one- 
dimensional Schrodinger equation, our method gives the well-known WKB J solution as the 
first adiabatic approximation. In the leading (second-) order approximation our approach 
leads to Bremmer's formula p!8[] for the over-barrier reflection (section 5). In section 6 
an expression for the transition probability due to a parametric excitation of a quantum 
oscillator is obtained from our results. Section 7 is devoted to a comparison of our 
calculations with exact solutions for a number of analytically solvable models. 



2 Dynamic equation on group manifolds 

We consider the linear operator (matrix) equation of the form 

dG/dt = B{t)G, G{to) = /, (7) 

where J is the unit operator and B(t) has a given time dependence. This is a pattern 
for a number of physical problems, including the one- dimensional Schrodinger equation 
and the spin precession in a time-dependent magnetic field. If the commutator algebra of 
operators B(t) is closed for all t, constituting a Lie algebra Q, one deals actually with the 
first-order dynamical equations on the Lie group G, generated by Q. Now the problem 
can be written in terms of the Cartan - Maurer one-form, 

dg g-^ = b{t)dt, git) G G, 6 G ^, (8) 

with the initial condition g{Q) = e - the unit element of G. Special problems are those 
where h{t) belongs to a Cartan subalgebra of i.e. h{t) C Q,Wt. In a case like that, 
Eq.(RD is integrated immediately. 



g{t) = exp 



to 



6(r)rfr 



e H, (9) 



where H is the corresponding Abelian subgroup of G. In general, (^ is a set of (non- 
linear) differential equations which cannot be reduced to quadratures. It is notable that 
the desired group element may be shifted by a properly chosen amount go{t), so the 
equation is rewritten in an equivalent form, g = go(t)gi 



Liiien 111 aii equivaieni luiiii, g = goyijgi, 
dgi = hi{t)dt, bi{t)dt = gQ^ (bdt - dgo g^^^^ go G Q. 



(10) 
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Thus the problem may be reduced to a more tractable one. 

Let us restrict ourselves to problems where b(t) approaches a Cartan subalgebra H 
asymptotically, as t — ±00, and evaluate the transition probability between eigen-states 
of operators representing Ti. The ^'-operator given by the following limit may be used 
(provided this hmit exists): 



S 



lim exp 

t,to—^^00 



B+{T)dT 



Gtoit) exp 



to 



111 



where B±{t) e H, and lim^ f^^-i-oo ll-B(t) — B±{t)\\ = 0. The probability of transition 
between the states given by the density operators P± a.t t ^ ±00 is 



W± = lim Tr 



'P+G{t)P_G\t)] = Tr (P+5P_5^) , 
since we assume that P± commute with B±. 



12) 



3 The adiabatic approximation 

At any given time t, the driving force b{t) G Q may be reduced to a Cartan subalgebra 
H, and the group element is decomposed respectively 

b{t) = v{t)(3{t)v{t)-\ (3{t) e n, (13) 
g{t) = w{t)h{t)w{ty\ h{t) e H. (14) 

Remarkably, if v has no t-dependence, that would fix the subalgebra Ti. for all t, and 
g would be obtained immediately, like in Eq.(^. We consider the problems where b{t) 
belongs to the Cartan subalgebra asymptotically, at t ^ ±00, so limf(t) = e. The 
equation resulting from (|^) would be 

w {dh h^^ + w~^dw — hw^^dwh^^^ = v(3v~^dt. (15) 

Splitting this equation to the subalgebras Ti and Q\H, we get a set of differential equations 
for h and w. We find an approximate solution of the Eq. ([T5|) for adiabatic processes, where 
the i-dependence of v is slow, the derivative dv/dt is small, and the condition 

\\v-'dv/dt\\ < WPW (16) 

is satisfied. Here the norm \\y\\ for an arbitrary element y of the Lie algebra Q is introduced, 

\\y\\ = v/Tr(rrt), yeg (17) 

and Y is the matrix belonging to the adjoint representation of the Lie algebra Q and cor- 
responding to the Lie algebra element y. When the condition (|I6|) holds, w is always close 
to V, and the deviation of g{oo) from the subgroup Ti is negligible. This is the meaning 
of the adiabatic theorem: the eigen-states of operators belonging to the subalgebra H are 
not subject to transitions. 
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The unknown group element may be replaced hj w = vexp{—uj), and it is assumed 
that u G G\'H. Small deviations from the adiabatic limit, producing non-adiabatic 
transitions, are obtained if we consider the first approximation in tu, which is expected to 
be of the order of v~^dv, discarding all higher-order terms. The result is 

dh h-^ - R{h){v-^dv - duj) = [13+ {uj(3 - (3uj)]dt, (18) 

where R{h)r] = hrjh'^ — f],^'r] ^ Q (note that R{h)ri = 0, if G 7Y). Separating the zero- 
and the first-order terms, we get two equations 



■ e H, (19) 



dh^ h^^ = (3{t)dt, so /iq = exp / (3{T)dT 

iJto 

R{ho)duj/dt + [hQ^dho/dt,uj] = R{ho)iv-^dv/dt). (20) 
The latter equation is also integrated immediately, 

R{hQ^)uj= f' R{hQ^){v-^v)dT, v = dv/dt\t=r- (21) 

Jto 

In the asymptotics, as soon as t, to ±oo, we get 



/oo 
R{hQ^){v-^v)dT, (22) 
-oo 



7 = 

and this element of Q determines the transition probability amplitude. In order to see 
that, let us insert the asymptotic value of the operator representing the group element 

g = V exp{—u!)h exp{u)v^^ ~ (e — u)ho{e + cu) ~ /io(e — R{hQ^)u) (23) 

in Eg . (|T^) for the transition probability. Assuming that P+ and represent different 
(orthogonal) eigen-states of operators corresponding to H, so that P+P- = = P-P+, 
one gets the following expression for the transition probability in the leading (second-) 
order 

W± = Tr(P+fP_f), (24) 

where F is the operator representing 7 G ^ in Eg. (^2]) . Note that any value may be taken 
for to in Eq-(19); changing it, say, to t'^, would result in a constant gauge substitution of 
7 for 7', 

rto 

(3{T)dT 



i = (K) S^o, K = exp 



(25) 



That would not change the probability in (|24D . The convergence of the integral in (|2^ ) 
depends on how fast the driving force b{t) is approaching its asymptotics in Ti. It is 
noteworthy that the present result extends the standard perturbation theory. If b{t) = 
(3q + A6i(t), where 61 (t) — > at ±00, then to the first order in A one has to set /io(T) = 
exp(/3o'r) in (0), and the result is an extension of the Born approximation. In general, 7 
indicates the deviation from the adiabatical limit. Even if the perturbation is not small 
absolutely, 7 may be small because of two different reasons: i) the change of v is slow, 
though it may be not close to unity, which is the case for small perturbations, ii) the 
deviation of h{t) from Ti takes place during a small time interval, and the integral is small 
result of that. 
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4 Spin-flip in a variable magnetic fleld 



The spin precession in a time-dependent magnetic field, the fundamental problem for 



NMR|T^, is determined by the Bloch equation for the spinor wave function, 

idi)/dt = {B ■ a)^, (26) 

where B = /iB(t), B(t) is a variable magnetic field vector, fi is the particle magnetic 
moment, and a are the Pauli matrices. The fundamental solution of Eq.(p6D is given by 
a unitary 2x2 matrix, so the group is SU{2). The matrix G is given by Eq.(0) with 
B = -i{B-a). 

Let us consider, for instance, the case where a pulse is applied in the x-direction, while 
the 2;-component is constant, B = {Bi{t), 0, -Bo}, and Bi{t) — as ^ ±oo. The adiabatic 
approximation holds if Bi{t) is a slow function of t. The elements which appear in ([T3| ) 
are 

Pit) = 2^[Bl + Bl{t)Y'^J,, v{t) = exp(0J2), (27) 

where tan 6 = —Bi(t)/BQ, and Ja is the basis in Q, represented by |(Ja. For this particular 
representation , we find from the Eq.(p!7D , that 



\\v~'dv/dt\\ = \dd/dt\- = ^i[Bl + Bl{t)Y/^ (28) 

and the general condition of applicability (0) leads to 

\de/dt\<^^l[Bl + Bl{t)Yl^ (29) 

The spin-flip, as given by Eq.(|2^), is determined by 7 = y4+_J2, and 

= r e^^^W^dr, a(r) = ^ [ [Bl + Bl{t)Y'Mt. (30) 

The spin probability is = |y4+_p. The accuracy of the approximation has been 

checked for a field where the exact analytical solution is available (see section 7). 

In the formula (^) dl{T) = Odr is the Riemannian length element on the path 
n(r) = (sin 6*, 0, cos 6). The geodesic curvature determined by Eq.(|), is equal to zero. In a 
more general case of the magnetic field configuration B = | i? | (sin ^ cos </>, sin 6' sin 0, cos^) 
the group element Qq leading to Eg . (p!0D will be chosen as Qq = exp(— 20J3). This trans- 
formation alters the phase Q;(r) in Eq.(|30|) to /^[sin^6' + (cos6' — ■^^)'^Y^'^\B\dt. An 

expansion of Q;(r) to the lowest-order non- vanishing in leads to Berry [g], Joye, Kunz 
and Pfister]^], Jaksic and Segert[|],|^ result for the spin-flip probability. 

5 Over-barrier reflection 

The Schrodinger equation, \E'" — f/(x)\E' = — A;^\E' is equivalent to the following first-order 
problem for the two-component function iplx), satisfying the equation 
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where 

Thus the coordinate x plays the role of the time parameter t. For the plane wave moving 
in the positive direction, = Cexp{ikx), so the upper component of ip vanishes, and 
reflection is like the spin flip. The problem of barrier penetration is represented by Eq.(|^, 
where 5^ = —aBa [a = is the diagonal Pauli matrix). Thus = aG^^a, the 
probability current j = —^ipaip = — |(\E''^' — \E'^'') is conserved, so we are dealing with 
the two-dimensional representation of the group G = SU (1,1). For U{x) decreasing 
rapidly ±oo, one has the transfer matrix 

I)' -If 1^ = 1. (32) 

The penetration probability amplitude is and the reflection probability is |6/ap. 

As soon as det B = k"^ — U{x) = P^ix), the Abelian subgroup is H = f/(l) in the region 
where p"^ > 0, and H = R under the barrier, where p"^ < 0. The 2x2 matrix diagonalizing 
B is 

y ^ f coshr] sinh?7 \ (2 ) = R f33) 

I sinhi] coshr/ J ' ^-^Pv ~ j^' 

Being applied to the plane wave , moving in the positive direction, the group element 
in the first-order adiabatic approximation g ^ vh^v'^ leads to the following expression 
for the wave function: 



^{x) = C 



(34) 



coshr]{x)e J^o'"^^' ^ + smhr]{x)e J^o'"^^' 

In the limit x +00 the parameter //(x) goes to zero, 

lim ^(x) = Cexp (i r pi^d^) (35) 

and the over-barrier reflection is absent in the first-order adiabatic approximation. Note 
that in the framework of this approximation the elements of the matrix V are considered 
to be slowly dependent on x, and exp(— r7(x)) = \J Thus a familiar WKBJ expression 
for the wave function: 

^{x) = Cn/^exp (t rpim) + C2\^exp {-i f p{m] . (36) 



P V Jxo 

is reconstructed. Remarkably, v belongs to a one-parameter subgroup of SU{1, 1), which 
makes the calculations simpler than in the general problem of spin precession (cf. section 
4). The over-barrier reflection is determined by (P^, where the element v is represented 
by the 2x2 matrix V , and -by the matrix: 



-i P pdx 

J^o=r I (37) 
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The probabihty of the over-barrier reflection is i? = |y4p, where 

1 2i r pdx u'{x) . 

A = - e J-o"^ ,„ V,; . dx 38 

The proposed method is valid when the unequality (|1^) is satisfied. This leads to the 
same condition of applicability as in the WKBJ approximation, 

\dp{x)/dx\ < p^{x) (39) 



Eq.(^) coinsides with the over-barrier reflection amplitude obtained by Bremmer|T! 
Bremmer's approximation was to divide a smooth potential into a large number of small 
layers. The momentum p{x) , being different in different layers, was assumed to be 
constant throughout a range of a particular layer. The ordinary WKBJ solution Eg . (|36|) 
was then obtained by discarding all reflections of the incident wave at any layer's boundary. 
Assuming that only single reflections at all boundaries of the layers take place, Bremmer 
found the over-barrier reflection amplitude. 

An advantage of our procedure is that leading to the same result (Eq.(^)) as Brem- 
mer's approach, our derivation does not demand any assumptions about qualitative char- 
acter of the wave reflection. Thus our method shows that the usual condition of the 
validity of WKBJ approximation (Eq.(^6D) is only needed in order to get Bremmer's 
formula. 

The reflection amplitude (^) may be compared with the one given by Maitra and 
Heller |T^. These authors use a perturbative approach with the WKBJ states as the 



unpertubed basis. According to Maitra and Heller the approximated reflection amplitude 
^M.H. is given by the matrix element of an effective potential between the usual WKBJ 
wave functions, i.e. 

^+oo g2j p p{y)dy 

Am.h.= Ueff{x,k) — dx (40) 

J-oo p{x) 

The effective potential of Maitra and Heller is given by the formula: 

n ( -3(pXx))^ , p"{x) 

Once Maitra and Heller use the perturbative arguments their expression for the reflection 
amplitude should be valid when the effective potential is small, i.e. 

\Ueffix,k)\<^k\ (42) 

Comparing the above conditions with that of applicability of our approximation we can 



see that our approach has a wider range of validity since inequality (^2]) follows from 
Eq.{^. 

In the limit of k"^ ^ U{x) the momentum p{x) becomes approximately equal to k. 
Neglecting U{x) in comparison with k"^ in the integral (|38[) , and integrating by parts, we 
obtain the result which corresponds to perturbation theory: 

1 r+oo 

Apert = ^ / ^'^^U{x)dx (43) 
ZlK, J — oo 
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6 Parametric excitation of a quantum oscillator. 



The parametric excitation of a quantum oscillator is the excitation of the oscillator under 
change of its parameters m = m{t) and Vt = Q{t). The general case with a time-dependent 
m(t) and may be easily reduced to m=const by changing variables t = J 

il' = mf2. The Schrodinger equation for the wave function of the quantum oscillator has 
the following form: 

.9-0 d'^ifj 1. 



dt 

For fl{t) the asymptotic conditions 

Q{t) 



t 



±oo 



are assumed. The assymptotic stationary states are: 



fx, fij 



2"n! V TT 



exp 



\HJ\ Vt+x 



(44) 
(45) 

(46) 
(47) 



The time-dependence of the quantum oscillator parameters m{t) and VL{t) allows for the 
transitions between different stationary states. 

A typical problem is to calculate the probability of transitions Wmn from the state 
il)n with the asymptotic 0^~^(x, t) at t — > — oo to the asymptotic state 0^^(x, t) at t — > 
+00. As is well known (e.g. Baz', Zeldovich and Perelomov||17||), in order to determine 
this probability of transitions Wmn it is sufficient to calculate the quantum mechanical 
coefficient 9 of the over-barrier reflection from the one-dimensional potential of a particle 
with momentum p(x) = ^Jk"^ — U {x) = fl{x), where fl{x) is the frequency function of the 
quantum oscillator. However, the analytic solution of the over-barrier reflection problem is 
known only for a number of special cases. When it is impossible to flnd an analytic solution 
for the problem of the over-barrier reflection, our approximated approach developed in 
section (^) may be applied. The approximated expression for the quntum mechanical 
coefficient 6 is given by the formula (B3), i.e.: 



2i n{t)dt Q,'{t)dt 



(4J 



When the parameter 9 is determined, the probability of transitions Wmn may be calculated 
using the Perelomov and Popov|]l6[ formula: 



Wmn 



n<! 
n>! 



m — n\ 



(49) 



where ?7,< 
tions. 



■min{m,n), n> = maxim, n) and P^{x) are the associated Legendre func- 
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7 Comparison with exact solutions. 
7.1 Spin precession in a magnetic field. 



The exact solution of the spin precession problem of Section 4 is known |[TT|1 for the mag- 
netic field 

B(t) = ;^f^^,0,/3o) (50) 
T \cosh(|;) ) 

Here [3q/T is the asymptotical precession frequency, T is the pulse duration. The appli- 
cability condition Eq. ( PD| ) allows for the application of our approach, when the inequality 

A//?o « 4^1 + ^0 (51) 



is satisfied. When /3o > 1, the inequality ( ^1|) holds for all /?i and the process would be 
adiabatic. Respectively, the perturbation theory can be applied for /3i -C /Sq- As known 
from the analytical solution of Eq. (|^ given in terms of the hyper-geometric function, 

= [sin(7r/5i)/cosh(7r/5o)]^. (52) 

Calculating the integrals in (|30|), it is usefull to change the variables as follows 

Pi/Po = tan/c, cos A; sinh(t/T) = sinh^. (53) 

Taking the integral for a with ri = 0, we have 

(^iO = Poi + Pi arctan(tan/ctanh^), (54) 

, f°° r /.NT tanhffi^ , , 

A+_ = sinfc/ sin[2a(0]-— j-^T%7TT (55) 
-'o (cosh^ - sm^ fc) 2 

The latter integral is reduced to a real form, as the pulse is symmetrical under the time 
inversion, so 9 is odd. The numerical calculation shows a wonderful accuracy of the 
approximation, namely, 

A+- ^ TT-TT 56 

cosh(7rpo) 

even for moderate values of (3q , in a wide range of f3i (see Figure 1). 

7.2 Over-barrier reflection for the potential U = j^^^- 

The analytic expression for the reflection amplitude is 



^ ^ sinh(7ra(l - VI -/j)) 

sinh(7ra(l + Vr^)) ^ ^ 

In the above formula the parameters a = k/'j and (3 = Uo/k"^ were introduced. Pertur- 
bation theory may be applied when k ^ Uq, i.e. < /3 <C 1. In that case the reflection 
probability is equal to 

P - 7T7- T 58 

^ 4sinh2 27ra ^ ' 
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Next, we consider the situation when fc^ > U{x). As follows from the inequality (|39|), in 
the cases when 

< 1 - /? (59) 

the over-barrier probability amplitude may be calculated by our method. It is suitable 
to change the variables z = e^^ and calculate the integral in the exponent of formula 
(p^). We obtain the following integral expression for the over-barrier reflection probability 
amplitude: 



(3 /•+°o 1 {2,/T^^{z + im-(3)z + l) + 2(1 - + 2 - /5) ^ dz 
(2^(z + l)((l-/3)z + l) + (2-/?)^ + 2) (i + z)((i-f3)z+l) 

A comparison of the exact and approximate (see Figure 2) probability amplitudes demon- 
strates a very good accuracy of the approximation (^), namely 

sinh7ra(l — ^/l + p) 



8 Conclusions 

In this paper we have proposed an adiabatic approach to the calculation of probabilities 
for quantum transitions. In the case when the one-parameter dependent Hamiltonian rep- 
resents a smooth curve in a Lie algebra, the original Schrodinger equation was interpreted 
as the dynamical equation on the corresponding group manifold. The main result of this 
work is expressed by Eq.(p2D that determines the Lie algebra element responsible for the 
non-adiabatic transitions. 

The problem of over-barrier reflection in one-dimensional quantum mechanics is very 
similar in our approach to the problem of spin-flip in a variable magnetic field (the differ- 
ence lies in the fact that for the over-barrier reflection problem the introduced evolution 
operator is an element of the group SU(l.l), and not of SU(2), as the spin evolution 
operator). 

We have tested our approach on simple problems for which approximate solutions are 
known. In the case of a spin in a time-dependent magnetic field our procedure leads 
to a spin-flip amplitude Eq . (|30|) . In the adiabatic limit Eq . (|30|) would coincide with 
the Berry]^], Joye, Kunz and PfisterQ], Jaksic and Segert[|],|0 result for the spin-flip 
probability amplitude. The application of our procedure to the over-barrier reflection gives 
Bremmer's formula[|r^ in the leading (second-) order approximation. It is remarkable 
that in order to obtain Bremmer's result, the usual condition of the validity of WKBJ 
approximation (Eq.(^)) is only needed. 

Being checked for two solvable models (the spin-flip in the Rosen-Zener magnetic 
field and the over-barrier reflection for the potential U = Uq/{1 + e'"^^) ), our adiabatic 
approximation not only gives the exponentially small character of the probabilities of 
the non-adiabatic processes, but completely describes the qualitative behavior of these 
probabilities as functions of the external parameters. The integrals (|S(]|), ( P^ show the 
same behavior under variation of the magnetic field amplitude (the amplitude of the 
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potential) as the exact solutions. It is interesting to note that in spite of the same 
condition of applicability as the WKBJ approximation, our approach is very successful in 
the calculation of the over-barrier reflection while the usual WKBJ approximation gives 
zero answer in all orders. The reason is that the WKBJ approximation is an asymptotic 
series that is unable to take the exponentially small variables into account. 
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